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Abstract. We present results of a monitoring campaign of the high-mass X-ray binary system 4U 1700-37/HD 153919, 
carried out with XMM-Newton in February 2001. The system was observed at four orbital phase intervals, covering 37% of one 
3.41-day orbit. The lightcurve includes strong flares, commonly observed in this source. We focus on three epochs in which the 
data are not affected by photon pile up: the eclipse, the eclipse egress and a low-flux interval in the lightcurve around orbital 
phase t/> ~ 0.25. The high-energy part of the continuum is modelled as a direct plus a scattered component, each represented 
by a power law with identical photon index (a ~ 1.4), but with different absorption columns. We show that during the low-flux 
interval the continuum is strongly reduced, probably due to a reduction of the accretion rate onto the compact object. A soft 
excess is detected in all spectra, consistent with either another continuum component originating in the outskirts of the system 
or a blend of emission lines. Many fluorescence emission lines from near-neutral species and discrete recombination lines from 
He- and H-like species are detected during eclipse and egress. The fluorescence Fe Ka line at 6.4 keV is very prominent; a 
second Ka line is detected at slightly higher energies (up to 6.7 keV) and a K/3 line at 7. 1 keV. In the low-flux interval the Fe Ka 
line at 6.4 keV is strongly (factor ~ 30) reduced in strength. In eclipse, the Fe Kfi/Ka ratio is consistent with a value of 0.13. In 
egress we initially measure a higher ratio, which can be explained by a shift in energy of the Fe K-edge to ~ 7.15 keV, which 
is consistent with moderately ionised iron, rather than neutral iron, as expected for the stellar wind medium. The detection 
of recombination lines during eclipse indicates the presence of an extended ionised region surrounding the compact object. 
The observed increase in strength of some emission lines corresponding to higher values of the ionisation parameter £ further 
substantiates this conclusion. 

Key words. Missions: XMM-Newton - Binaries: massive, eclipsing - Stars: individual: 4U1700-37/HD153919 - accretion - 
scattering 



1. Introduction 

4U 1700-37/HD 153919 is a high- mass X-ray binary (HMXB) 
discovered by the Uhuru satellite (Lfones et al.lll973l) . HMXBs 
are relatively young (< 10 7 year) systems composed of a mas- 
sive OB-type star and a compact object, either a neutron star 
or a black hole (for a recent catalogue of X-ray binaries see 
Liu et al . 2000). X-rays are produced when matter from the OB 
star accretes onto the compact object. The majority (~ 80%) of 
the HMXBs are Be/X-ray binaries. These systems are often X- 
ray transients, only producing X-rays when the compact object 
accretes matter by passing through the dense equatorial disc of 
a Be-star. The other group o f systems hos t a more massive OB 
supergiant (for a review see lKaperll200ll) . About a dozen OB 
supergiant X-ray binary systems are known in the Milky Way 
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and the Magellanic Clouds. Among them are a few that host a 
black-hole candidate (e.g. Cyg X-l). 

Accretion takes place through either the strong stellar wind 
of the optical companion or Roche-lobe overflow. In a wind- 
fed system the X-rays ionise a relatively small region in the 
stellar wind due to the low X-ray luminosity (Lx ~ 10 35 - 10 36 
erg s _1 ). The size of the ionisation zone can be determined 
from the orbital modulatio n of UV resonance lines forme d 
in the wind JHatchett & McCravl[l977l Ivan Loon et alJl200 lh . 
In a Roche-lobe overflow system, matter flows via the inner 
Lagrangian point to an accretion disc. A powerful X-ray source 
(Lx ~ 10 38 erg s~') is produced, which ionises almost the 
whole stellar wind, such that only in the X-ray shadow of the 
optical companion a stellar wind can emerge (e.g. LMC X-4; 
iBoroson et alJl99^). For a more g en eral review on X-ray bina- 
ries we refer to e.g. lNagasd ( ll989l) or lLewin et alJ dl995h . 
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X-ray spectra of HMXBs are often described by a power 
law with photon index a ~ 1, modified at high energies (above 
10-20 keV) by an exponential cutoff. A spectrum of this form 
can be produced by inverse Compton scattering of soft X-rays 
by hot el ectrons in the ac cretion column near the compact ob- 
ject (e.g. Hav akawall985h . This continuum component is often 
scattered by the stellar wind of the optical companion. This re- 
sults in a similar component, but with a different absorption 
column depending on the orbital phase of the system. In the 
high energy part (> 10 keV) of the spectrum of most X-ray 
pulsars broad absorption lines have been discovered, which are 
attributed to cyclotron resonant scattering. From the energy of 
these lines strong magnetic fields are inferred (B ~ 10 12 Gauss; 
e.g. Trump er" ;t alll978l) . 

In some HMXBs a soft excess at ~ 0.1 — 1 keV is detected 
of which the physical origin is not well understood. If this com- 
ponent has a thermal origin, its relatively low temperature re- 
quires in some cases a size of the emission region much larger 
than the size of the compact object, of the order of the mag- 
netospheric radius, where the pressure of the infalling matter 
balances the pressur e of the magnetic fie l d of the compact ob- 
ject (e.g. Cen X-3: iBurderi et alJl200ol) . Iffickox etafl J2004) 
show that a soft excess is very common in X-ray binary pulsars 
and that observations that do not show evidence of a soft excess 
are probably li mited by high a bsorption, low flux or insufficient 
soft sensitivity. Hick ox et alJ (J2004) find that the soft excess in 
luminous (Lx > 10 38 erg s _1 ) sources can be explained by re- 
pocessing of hard X-rays from the neutron star by optically 
thick, accreting material. For less luminous (Lx < 10 36 erg s _1 ) 
sources the soft excess is probably due to other processes, e.g. 
emission from photoionised or collisionally heated diffuse gas 
or thermal emission from the surface of the neutron star. 

With the advent of the X-ray observatories Chandra and 
XMM-Newton it has become possible to study the X-ray spec- 
tra of HMXBs in unprecedented detail. Their eclipse spectra 
show many emissio n lines and radiative recombination con - 
tinua (e.g. Vela X-l lSchulz et al.l2002UGoldstein et all2004 . 
Such an emission line spectrum originates in the ionisation 
zone surrounding the X-ray source and becomes detectable 
during eclipse due to the supression of th e strong con t inuum 
emission produced by the X-ray pulsar. IPaerels et alJ (^000) 
show that the Chandra spectrum of Cyg X-3 is consistent 
with a cool, optically thin, photoionised plasma in equilib- 
rium, whereas the spectrum of 4U 1700-37 shows a more "hy- 
brid" plasma in which also evidence is found for a collis ionally 
ionised plasma (B oroson et alJl2003l) .[Woidows ki et alJ (J2003) 
report on an apparent change of plasma properties from eclipse 
to eclipse egress of Cen X-3, which they can explain with a 
photoionised plasma for which the e ffects of resonant line scat- 
tering are taken into account as well.fwatanabe et al. ( 2003J re- 
port on a fully resolved Compton-scattered Fe Ka fluorescence 
line in GX 301-2. They show that the flux ratio of the line and 
shoulder together with the energy distribution of the shoulder 
can be used to directly constrain the absorption column, tem- 
perature and abundances of the scattering medium. 



1.1. 4U 1700-37 

Located at a distance of 1.9 kpc jAnkav et alJ l200lh 
4U 1700-37 i s powered by the dense stellar wind (M~ 
10~ 5 M s yr" 1 ; IClark et alJ 120021) of the Q6.5 Iaf+ super- 
giant HD 153919 Jjones et alJll973t iMason et alJll97fih . It is 
the hotte st and most l umino us optical companion known in 
HMXBs. lAnkav et alJ i l200ll) propose that the system escaped 
the OB association Sco OBI ~ 2 million years ago, travelling 
with a velocity of ~ 75 km s _1 . They show that the progeni- 
tor of the compact object should have evolved into a supernova 
within ~ 6 million years (if it was born as a member of the as- 
sociation). T his corresponds to a lower limit of its initial mass 
of ~ 30 M© iSchalleret alJl992l) . 

Period ic X-ray pulsations h ave not been detected so far 
(although iBoroson et al.l [2003 report quasi-periodic oscilla- 
tions), so_JhatJhe_Mtare_of the compac t object remains un- 
clear JGottwald et all 1 19861 IClarketalJ 12002). Most likely, 
4U 1700-37 is a neutron star since its X-ray spectrum 
is well desc ribed by the standard accreting pulsar model 
dWhiteetalJ fT983). There are indications that the mass of 
the compact obje ct is larger than 2 Mq (Ru bin et alJll99rj: 
IClark et all l2002h. If true, this m ight qualify it as a low- 
mass black-hole ( Brown et al. 1996) or a massive neutron star, 
like in Vela X-l JBarziv et alJl200ll lOuaintrell et~ai1l2003l) . 
iRevnolds et alJ i 19991) report the possible presence of a cy- 
clotron feature at ~ 37 keV in BeppoSAX observations. If real, 
4U 1700-37 has a magnetic field B ~ 2.3 x 10 12 Gauss, char- 
acteristic of a HMXB accreting pulsar. 

This system is unique in showing broad and v ariable emis- 
sion li nes in the ultraviolet part of the spectrum. iKaper et"all 
(1990) propose that these lines are produced by photons from 
strong EUV emission lines close to the He n Lyman-/? transition 
(256 A), which are Raman scattered by He n ions in the stellar 
wind of the optical companion via the corresponding down- 
wards transition near the Hen Balmer-a line at 1640 A. The 
input EUV emission lines can not be studied direct ly due to the 
opaqu eness of the interstellar medium. However, IKaper etafl 
(1990) show that these emission lines originate close to the 
compact object, causing the orbital modulation of the Raman 
scattered lines. Although the central wavelength of the input 
EUV emission lines is known, they could not be identified due 
to the multitude of candidate line features in this part of the 
spectrum. 

The X-ray spectrum of 4U 1700-37 can be modelled as 
a power law with a high-energy cutoff JHaberl et alJ Tl989: 
IRevnolds et al J 19991) . modified by an absorption column which 
depends on the orbital phase of the sy stem. A sharp in- 
crease in absorptio n at late orbital phases (Ma son et al.lll976l 
iHaherl etal.lll 989) is explained by the presence of a region of 
enhanced density trailing the compact obje ct in its orbit, such 
as a photo-ionisation wake (Blondin et al. 1990). The presence 
of such a wak e in 4U 1700-37 h as been confirmed with optical 
spectroscopy (Kap er et all 1994b . A similar increase in absorp- 
tion has b een reported for other HMXBs, like Vela X-l and 
Cen X-3 faaberl & White! 1990tlNagase et alll992l) . 

A soft excess was found in Ginga observations 
JHaberl fe Davl 1 19921) and modelled as a bremsstrahlung 



A. van der Meer et al.: XMM-Newton X-ray spectroscopy of the high-mass X-ray binary 4U1700-37 at low flux 



3 



Phase 0.22 - 0.31 



51960.7 




51957.8 51957.9 51958.0 51958.1 

MJD (days) 



51958.2 




51958.9 51959.0 
MJD (days) 




51959.6 51959.7 51959.8 51959.! 

MJD (days) 







Phase 0.07 - 0.17 




eclipse 


egress 









51960.8 



51960.9 51961.0 
MJD (days) 



51961.1 



Fig. 1. 0.3-1 1.5 keV lightcurve of 4U 1700-37 obtained with 
the EPIC/MOS-2 camera at the end of February 2001 in four 
orbital phase intervals. The data have been extracted from a 
circular region centered on the source and binned in time-bins 
of 10 seconds. The data have not been corrected for pile-up, 
which mainly distorts the shape of the spectrum but may also 
affect the value of the high count rates, e.g. during flares. We 
selected three intervals in the lightcurve (indicated by the ver- 
tical dashed lines), in which the spectrum does not suffer from 
pile-up, i.e. the eclipse, the eclipse egress and a low-flux inter- 
val, centered on orbital phase ~ 0.25. 



component with a temperat ure of kT ~ 0.5 ke V. In a second 
set of Ginga observations Haberletal ] (11994 also found a 
difference in temperature before and after eclipse, which they 
tentatively explain with a bow shock in front of the compact 
object, heating the stellar wind to a higher temperature. 



Table 1. Log of XMM-Newton observations of 4U 1700-37 
in the period 17-20 February 2001 . To calculate the or- 
bital phase we used the ephemeris of Rubi n et al or- 
bital period P 01 -b = 3.411581(27) d and mid-eclipse To = JD 
2448900.873(2). 
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iBoroson et alJ J2003I) obtained Chandra observations of 
4U 1700-37 in August 2000 in the orbital phase interval 
0.65 < <p < 0.8. They find emission lines in the 4-13 A 
range, which are identified with fluorescence emission lines 
from near-neutral species and recombination lines from He- 
and H-like species. The X-ray source was observed during an 
intermittent flaring state, in which the strengths of the lines 
vary. The triplet structure of Si xra and Mg xi is resolved but 
weak, but could be used to tentatively identify the plasma as 
a low-density hybrid of a photo-io nised and a collisionally 
ionised plasma ( IBoroson et alJ2003l) . 

We monitored 4U 1700-37 with XMM-JVewfon around the 
four quadratures in the orbit, including part of the eclipse and 
the eclipse egress. The goal of our observing campaign is to 
study the X-ray spectrum as a function of orbital phase and 
to derive information on the physical origin of the continuum 
and line emission. In Sect. 13 we describe the observations and 
the techniques that we used to reduce the dataset; in Sect. [3] 
we present the spectral analysis, and in Sect. 0]we discuss our 
results. 

2. Observations 

We have obtained fo ur 20-30 ks observations of 4U 1700-37 
with XMM-Newton djansen et alJ l'2001) in the period 2001, 
February 17-20. The log of observations is listed in Tabled 
Fig. [2 shows the l ightcurve of the four intervals obtained 
with EPIC/MOS-2 dTurner et alJl200lh. To determine th e or- 
bital phase we used the ephemeris of Rubi n et alJ ( 1996). The 
X-ray flux of 4U1700-37 is strongly variable (factor ~ 10). 
Numerous flares are observed with a typical duration of the or- 
der of half an hour. 

Due to the brightness of the optical companion (My = 
6.5), we were forced to switch off th e Optical Monitor . 
The EPIC/MOS-1 and EPIC/PN detectors ( Struder e t alJ200lb 
were used in timing and burst mode, respectively, to cope with 
the high count rate of the X-ray source. EPIC/ MOS-2 was used 
in sm all-window mode. For the RGS gratings Jden Herder et alJ 
l200lh default settings were used. 

We reduced the data using the Science Analysis System 
(SAS) version 5.3.3. Since this version does not yet con- 
tain a developed task to reduce EPIC/MOS-1 timing and 
EPIC/PN burst mode, in this paper we will concentrate on the 
EPIC/MOS-2 and RGS datasets. The results from the EPIC/PN 
and EPIC/MOS-1 will be reported in a forthcoming paper. 
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Fig. 2. XMM/EPIC observation of 4U 1700-37 during eclipse egress (0.14 < <p < 0.17). The continuum emission is modelled 
as three power laws, representing a direct (solid line), a scattered (dashed-dotted line) and a soft component (dotted line). The 
spectrum clearly shows several fluorescence emission lines and some discrete recombination lines of He- and H-like species (see 
Table[3]for their identification). Note that the data in this plot are folded with the response matrix of the instrument. Therefore, 
some of the features (e.g. at 1.8 keV) in the plot are of instrumental origin. The lower panel shows the residuals of the model to 
the spectrum. 



Since the pipeline products were created when the reduc- 
tion software was still in an early stage, we created new event 
lists for the EPIC/MOS-2 data. These event lists were filtered 
following the data analysis threads provided by the XMM- 
Newton team. Because the instrument was set in small-window 
mode (to reduce the read-out time) our strong source produced 
numerous events over almost the whole frame, which prevented 
us from measuring the background. Therefore, we extracted 
a background spectrum using special event lists provided by 
the XMM-7Yewfo« team. 1 These are long exposure observa- 
tions of fields with a typical background and without any bright 
sources. Since our source is bright, the background contribu- 
tion is negligible (~ 0.1% in the energy range 0.3-11.5 keV 
and < 2% in the energy range 0.3-2.0 keV). 

Due to the high count rate a large fraction of the dataset is 
affected by pile-up. The basic concept of pile-up is that a sat- 
uration effect takes place in which two or more separate pho- 
tons fall onto the same pixel within a read-out cycle and are 



1 http://xrnrn.vilspa.esa.es/external/xrnrn_sw_cal 



counted as one with a total energy equal t o the sum o f their 
energies (for a detailed description see e.g. lBalletlll999l) . This 
effect distorts the source spectrum and the measured flux. The 
SAS software provides the tool EPATPLOT to analyse the dis- 
tribution of patterns produced on the CCDs by the incoming 
photons. Two or more photons counted as one have a lower 
probability to produce a single-pixel pattern than one photon 
with the same energy. This affects the observed pattern distri- 
bution, which can be compared with the expected distribution 
to assess the pile-up. However, with this method it is not clear 
to what extent the spectral parameters (derived from fits to the 
extracted spectrum) are affected by pile-up. 

In order to assess the magnitude of the distortion and to 
eliminate the part of the point spread function (PSF) affected 
by pile-up, we created multiple event lists, extracted from an- 
nuli centered on the source, with increasing inner radii (with 
steps of 2.5") up to a maximum of 25". We modelled the spec- 
tra with the same basic model of two absorbed power laws (see 
Sect. 13. II . The best-fitted spectral parameters were then ex- 
amined as a function of inner radius of the extraction region. 
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One expects the parameters to stabilise beyond the radius at 
which pile-up no longer affects the remaining part of the PSF. 
However, we find that the parameters continue changing, some- 
times even up to the maximum inner radius, while EPATPLOT 
already gives satisfactory results. We conclude that the model 
of the energy dependence of the PSF is not known accurately 
enough for our purpose. 

In the current study we concentrate on data with a suffi- 
ciently low countrate. We select three intervals in the lightcurve 
(see Fig.Q: eclipse, eclipse egress and a low-flux part centered 
on orbital phase d> ~ 0.25. These intervals of low flux are not 
uncommon in 4U 1700-37 (e.g. IBoroson et alj|2003l) and are 
also observed in other HMXBs, e.g. 4U 1907+097 studied by 
lin 'tZand et al 1 dl997l) . who refer to this behaviour as "dipping 
activity". 

We also analysed the RGS data extracted in the same time- 
frames. However, the source is highly absorbed, resulting in 
a very low RGS-countrate. After rebinning the data to 20-25 
counts per bin the RGS spectra were used to verify the best-fit 
models for the EPIC/MOS-2 data. 

3. Spectral Analysis 

The EPIC/MOS-2 camera covers the energy range 0.3- 
11.5 keV. In previous studies based on a w ider energy range 
(EXOSA T. 3-20 keV. IHaberl et al.1 fl989l BeppoSAX, 0.5- 
200 keV. lRevnolds et al.ll 19991) the continuum of 4U 1700-37 
is modelled as an absorbed power law with photon index a, 
modified above a break energy E c by a high-energy cutoff of 
the form exp KEV - E )/Ef]. 

ISako et al.( <fl999) obtain a good fit of the ASCA eclipse 
spectrum of Vela X-l using two absorbed power laws with 
identical photon index, but different column density Nh and 
normalisation to model the continuum. One power law de- 
scribes the direct component, that originates from close to the 
neutron star. A part of this component is scattered by the ex- 
tended stellar wind of the optical companion and is described 
by the second powerlaw. lBoroson et al apply the same 

method to fit the Chandra spectrum of 4U 1700-37. We follow 
the same strategy modelling the XMM-Newton spectra. We ap- 
ply three models to fit the continuum at eclipse, egress and low 
flux. The emission lines are subsequently fit by gaussians. The 
spectra have been analysed using XSPEC version 1 1.2.0, after 
rebinning the spectra to 25 counts per bin. 

3.1. Continuum Emission 

We first describe the eclipse egress spectrum, where a clear dis- 
tinction can be made between the different continuum compo- 
nents (see Fig. |2). We discuss the high-energy and the low- 
energy part of the spectra separately. Note that the mentioned 
parameters correspond to the models fitted to the full 0.3-1 1.5 
keV energy range. Table[2]lists the models we used to fit these 
three spectra of which the parameters are constrained: model A, 
which consists of three power laws, model B with two (or three) 
power laws modified by a high-energy cutoff, and model C, 
similar to model A, but with a blackbody component replacing 
the power law describing the soft component. 



3.1 .1 . Direct and scattered component 

We model the high energy part of the eclipse egress spectrum 
as two power l aws, each of them modified by different absorb- 
ing columns (Morrison & McCammon 1983) depending on the 
orbital phase of the system (for the egress spectrum these are 
N H ~ 8.3 x 10 23 cnT 2 and Nh ~ 1.8 X 10 23 cirT 2 ; see model A 
in Table |2j. Both components have the same power law in- 
dex a ~ 1.4. The hard, direct component is assumed to orig- 
inate from the near surroundings of the compact object and 
the scattered component is produced by Thomson scattering 
by electrons in the extended stellar wind of the optical com- 
panion. Since Thomson scattering is energy independent, the 
photon index of the scattered component is expected to be the 
same, and is therefore fixed to the index of the direct compo- 
nent. The best-fit values of the model parameters are given in 
Table |2 Both components are also detected during the part of 
the eclipse covered by our observations (0.07 < <p < 0.14), al- 
though the flux of the direct component is reduced by a factor 
~ 45 compared to egress. 

In the low-flux spectrum it is difficult to separate the two 
components. The model produces comparable strengths and 
absorption coefficients (Model A in Tabled} with a photon in- 
dex that is slightly higher (a ~ 1.9) than in eclipse and egress. 
Modelling the spectrum with only one power law modified by 
a high-energy cutoff gives an equally good fit (model B in 
Table |3- The resulting photon index is lower (a ~ 1.1) than 
in model A. The cutoff energy in this model is ~ 5.4 keV 
(see Table Et, which is in ag reement with previous observa- 
tions (e.g. IHaberl etalJll989l who find E r = 6.6 ± 0.7 keV 
in EXOSAT observations and Reynolds et alJll999l who find 
E c — 5.9 + 0.2 keV in BeppoSAX observation s). However, 
the folding energy we find (E f ~ 1 1 keV) is low; IHaberl et all 
dl989l) find E f = 21 + 4 keV and lRevnolds et alJ dl999l) mea- 
sureE/ = 23.9 + 0.9 keV. 

The problem with applying this model to our dataset is that 
the folding energy is outside the energy range of EPIC/MOS-2 
(0.3-11.5 keV) and its value is easily influenced by the value 
of the photon index. If we apply the model to the eclipse and 
egress spectrum, a high-energy cutoff yields an unacceptable 
fit, and the fit parameters can not be constrained. If we fix E c 
and E f to the values obtained from the low-flux spectrum, the 
fit parameters can be constrained (see Table|2j. This is also the 
case if w e fix their values to th e values of the BeppoSAX obser- 
vations dRevnolds et al.ll999l) . The fits with model A and with 
this "fixed" version of model B are of comparable quality for 
the eclipse spectrum, but for egress fits with model B are worse 
than fits with model A (x 2 /d.o.f. = 578/522 for model B com- 
pared to 558/522 for model A; d.o.f. = degrees of freedom). 
Note that model B also indicates a higher photon index in the 
low-flux interval compared to eclipse and egress. 

3.1.2. Soft Excess 

After having obtained the best fit for the models discussed 
above, an excess remains in all cases in the low energy part 
of the spectrum (see Fig. [3}- A soft excess is reported for many 
X-ray spectra of HMXBs, but the nature of this component is 



Table 2. Parameters used to describe the continuum models of the eclipse egress, eclipse and low-flux spectra. All errors and upper limits are a 1 cr confidence level. Three 
different models are considered. Model A consists of three power-law components with identical photon index, but different normalisation, which are modified by different 
absorption columns. Model B contains two power-law components, modified by a high-energy cutoff. Model C includes a blackbody component replacing one of the three 
power-law components of model A. 





egress 


eclipse 


low-flux 




Model A Model B Model C 


Model A Model B Model C 


Model A Model B Model C 


direct component 
N H (xlO 22 cm- 2 ) 

Normalisation (xl(T 2 photons s -1 cirr 2 keV -1 at 1 keV) 
Unabsorbed flux (xlCT 10 erg s -1 cnT 2 ) 


83 ± 2 85 + 3 83 +2 
44 ± 6 12 ± 2 44 ± 6 
49 ± 3 34 ± 5 49 ± 3 


59 + 8 89+ 12 61 +9 
1.1 +0.4 0.47 + 0.19 1.1 
1.1+0.3 1.0 + 0.4 1.2 + 0.4 


20 + 3 6.8 ± 0.4 17 + 2 
8.7 + 2.4 3.8 + 0.6 9.3 + 2.0 
5.3+^ 5.3 ± 1.0 6.1 ± 1.0 


scattered component 
N H (xlO 22 cm" 2 ) 

Normafisation (xf(T 2 photons s -1 cirr 2 keV -1 at 1 keV) 
Unabsorbed flux (xKT 10 erg s~' cnT 2 ) 


17.5 ±1.3 18.7 + 1.6 17.6 + 2.2 

1.6 + 0.3 0.67 + 0.13 1.6 + 0.4 

1.7 + 0.2 1.9 + 0.4 1.8 + 0.3 


11.6 + 0.7 11.2 + 0.7 12.3 + 1.0 
0.56 + 0.12 0.24 + 0.05 0.61+0.15 
0.60 + 0.06 0.54 + 0.12 0.66 + 0.07 


7.2 + 0.3 ... 6.4 + 0.4 
8.2+1.4 ... 6.3 + 1.4 
5.3+0.7 ... 4.1+0.7 


soft component 
N H (XlO 22 cm" 2 ) 

Normalisation (xl0~ 4 photons s -1 cirr 2 keV -1 at 1 keV) 
Blackbody temperature (keV) 
Normalisation (xlO 32 erg s -1 kpc~ 2 ) 
Unabsorbed flux (xlO"' 2 erg s -1 cirr 2 ) 


0.56 + 0.09 0.36 + 0.08 0.30 + 0.17 
4.3+0.5 2.6 + 0.3 

74+ - 34 
2 1 +1J 

4.7 + 0.5 7.5+0.8 1.8^ 


0.24 + 0.04 0.13 + 0.04 < 0.05 
2.3 + 0.2 1.7 + 0.2 

77+°- 11 

"•"-0.03 

1.3+0.4 

2.5 + 0.3 3.7 + 0.4 1.1+0.3 


0.43 + 0.08 0.32 + 0.07 1.0 + 0.4 
4.7 + 0.9 3.4 + 0.7 

0.13 + 0.04 
< 270 

3.0 + 0.5 4.7+1.1 <230 


Cutoff energy (keV) 
Folding energy (keV) 

Photon index 

^ 2 /d.o.f 


5.4 fixed 
10.6 fixed 

1.35 + 0.05 0.61+0.05 1.35 + 0.05 
558/522 578/522 559/521 


5.4 fixed 
10.6 fixed 

1.36 + 0.12 0.77 + 0.12 1.37 + 0.13 
444/408 445/408 443/407 


5.4 + 0.4 

1.87!™* 1.08 + 0.09 1.83 + 0.08 
485/390 487/390 475/389 
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Fig. 3. XMM/EPIC observation of the 4U 1700-37 during the low-flux interval. A soft excess below ~ 2 keV is clearly present. 
The solid line represents a model that consists of one absorbed power-law component. The lower panel shows the residuals of 
the model to the spectrum. 



not clear. To investigate the physical origin of this component 
we fit the soft component with four different models, which are 
discussed below: bremsstrahlung, blackbody, electron scatter- 
ing and dust scattering. In all cases it is found that the model 
component needs to be modified by a moderate absorption col- 
umn, close to the interstellar value Nh ~ 0.3 x 10 22 cm" 2 for 
this source, based on th e observed reddening E(B - V) = 0.59 
fautchingset all 1973b and N H /E(B - V) = 0.58 x 10 22 cirT 2 
mag" 1 llKilianll992l) . 

Adding a brems strahl ung component, a s proposed by 
lHaberl & Davl dl992h and lHaberl et all dl994l) . improves the 
fit to all three spectra significantly (^ 2 /d.o.f. = 476/389 
for the low-flux spectrum). However, the temperature of the 
bremsstrahlung component can not be constrained in any of 
the spectra. Actually, this bremsstrahlung component tends 
to mimic the shape of a power law. This suggests that 
bremsstrahlung is not the physical origin of the soft excess. 

Fitting this excess with a blackbody component gives a 
similar improvement of the fit (j^ 2 /d.o.f. = 475/389 for the 
low-flux spectrum). For this model the temperature kT ~ 
0.13 keV. The flux is not constrained in the low-flux spectrum, 



-2 



but we can place an upper limit of 2.3 x 10 10 erg s 1 cm 
In the eclipse and egress spectra the resulting temperature is 



slightly higher (kT ~ 0.77 keV and kT ~ 0.74 keV, respec- 
tively). The unabsorbed flux is constrained to 1 . 1 x 10" 12 erg s" 1 
cm" 2 and 1.8 x 10~ 12 erg s" 1 cm" 2 , respectively (see model C 
in Tables- 
Modelling the soft excess with yet another (third) power 
law with a photon index fixed to that of the direct and 
scattered component also gives a significant improvement 
(X 2 /d.o.f. = 485/390) and all the parameters are constrained 
(see Table |5Jl. This component could, like the scattered com- 
ponent of model A, represent a Thomson scattered version of 
the hard, direct component, with a very low absorption column 
(N H ~ 0.43 x 10 22 cm" 2 ). 

Scattering by dust grains becomes a possibility at larger dis- 
tance from the system (where the absorption is low) as found by 
Ebis awa et alJ dl996l) in the ASCA observations of Cen X-3. 
We add 2 to the photon index of the soft component compared 
to the index of the hard component, since dust scattering, un- 
like Thomson scattering, is energy dependent (oc I IE 2 ). This 
model results in an equally good fit as Thomson scattering for 
the low-flux and egress spectrum. Since the difference between 
Thomson scattering and dust scattering in our models is pro- 
duced by a change in the photon index, this will only be visible 
as a difference in flux at the high energy part of the spectrum. 
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Fig. 4. XMM/EPIC observation of 4U 1700-37 during eclipse. The spectrum clearly shows several fluorescence emission lines 
and some discrete recombination lines of He- and H-like species with their identifications. The continuum model used for this 
spectrum is model A (see Table |2}; the RGS spectrum is also included. The EPIC/MOS-2 residuals are shown in the middle 
panel, and the RGS residuals in the lower panel. 



However, this part of the spectrum is dominated by the flux 
of the direct component and therefore the difference between 
dust or Thomson scattering is difficult to measure. For the 
eclipse spectrum, where the direct component does not dom- 
inate the high-energy part of the spectrum, the dust scattering 
model yields a significantly worse fit Or 2 /d.o.f. = 453/408) 
than Thomson scattering (x 2 /d.o.f. = 444/408). 



3.2. Line Emission 



Following the approach of Bo roson et al.N2003) we also try 
to model the soft excess as a blend of gaussians only. We delete 
the soft continuum component from model A and refit the data. 
At low energies the fit now contains broad lines with a width on 
the order of 150-180 eV, whereas the resolution of EPIC/MOS- 
2 has a width of ~ 40 eV at 1 keV. The absorption of the direct 
and scattered components is "lower" in eclipse than in egress 
(N H ~ 5.6 x 10 22 crrT 2 and N H ~ 30 x 10 22 crrT 2 in eclipse 
w.r.t N H ~ 15.3 x 10 22 cirT 2 and N H ~ 82 x 10 22 cirr 2 ). 

We also investigated the soft component in the RGS dataset 
(see Hig - . Due to the low RGS count rate it was not possible 
to set further constraints. 



As in other HMXBs (e.g. Vela X-l and Cen X-3: ISako etail 
1999; Ebisawa et al. 1996), 4U 1700-37 shows many emission 
lines in the eclipse and egress spectrum, when the continuum 
emission is at a minimum. Only some of these lines can be 
detected in the low-flux spectrum. We fit these lines as gaus- 
sians and list them in Table[3]if their significance is larger than 
2cr. For the lines detected in some, but not all spectra, we set 
an upper limit on the line flux in case of a non-detection. We 
do not correct for absorption in calculating the flux, since the 
geometry of the line forming region is not known. 

As the energy resolution of the EPIC/MOS-2 instrument is 
not sufficient to resolve all lines, many of them are blended. 
This complica tes the fine identification process (for which we 
use Jo hnson & Sofa ll985: Drake 1988), especially in the low 
energy range. The lines that we can identify are fluorescence 
emission lines from near-neutral species and discrete recombi- 
nation lines from He- and H-like species. We will concentrate 
on the iron complex and the lines that are clearly resolved and 
report on their behaviour, moving from eclipse into egress and 
their presence in the low-flux spectrum. 
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Table 3. List of the emission lines that have been detected in the eclipse, the eclipse-egress and the low flux spectra. All errors 
and upper limits are 1 cr. If a line is detected in one spectrum but not in another, an upper limit on the flux is given. 



eclipse 



egress 



low-flux 



eclipse 



egress 



low-flux 



Energy (keV) 
Width (eV) 

Flux (10" 5 photons cnT 2 s" 1 ) 
Identification 

Energy (keV) 
Width (eV) 

Flux (10" s photons cnT 2 s" 1 ) 
Identification 

Energy (keV) 
Width (eV) 

Flux (10" 5 photons cm" 2 s" 1 ) 
Identification 

Energy (keV) 
Width (eV) 

Flux (10" 5 photons cm" 2 s" 1 ) 
Identification 

Energy (keV) 
Width (eV) 

Flux (10" 5 photons cm" 2 s"') 
Identification 



Energy (keV) 
Width (eV) 

Flux (10" 5 photons cm" 2 s"') 
Identification 

Energy (keV) 
Width (eV) 

Flux (10" 5 photons cm" 2 s"') 
Identification 

Energy (keV) 
Width (eV) 

Flux (10" 5 photons cm" 2 s" 1 ) 
Identification 



0.86 ± 0.02 0.87 ± 0.02 
71 ± 19 84 + 32 

4.0+1.0 5.5 ±2.0 
NeKa, OvmRRC 

1.01 ±0.01 1.01 fixed 

< 13 13 fixed 

1.0 ±0.3 <0.52 
Ne x Lya 

1.30 ±0.01 1.32 ±0.03 
77 ± 17 < 87 

3.1 ±0.7 2.2 ±0.7 
Mg Ka , Mg XI Hea 

1.47 ±0.02 1.47 fixed 

< 31 31 fixed 
0.6 ±0.2 <0.78 
Al Ka , Mg xil Lya 

1.75 ±0.01 1.77 ±0.01 

< 26 < 26 
5.0 ±0.6 4.8 ±0.9 
Si Ka , Al xiii Lya 

1.85 ±0.01 1.86 ±0.06 

< 20 < 101 
3.4 ±0.8 1.6^j 
Si xm Hea 

2.00 ±0.01 2.02 ±0.01 

< 18 <31 

1.6 ±0.3 1.9 ±0.5 
Si xiv Lya 

2.33 ±0.01 2.32 ±0.01 

34 ± 11 < 34 

6.3 ±0.7 6.1 ±1.4 
S Ka, AlxmRRC 



0.86 ± 0.02 
60 ±27 
5.3 ± 1.8 

1.01 fixed 
13 fixed 

< 1.7 

1.31 fixed 
87 fixed 

< 3.3 

1.47 fixed 
31 fixed 

< 0.15 

1.76 fixed 
26 fixed 

< 0.35 

1.83 ±0.03 

< 42 
2.7 ± 1.5 

2.01 fixed 
20 fixed 

< 4.2 

2.33 fixed 
34 fixed 

< 4.9 



0.01 



2.46 
< 18 
1.6 ±0.4 
S xv Hea 



9 47+0.04 

'"^ -0.02 

115 ±29 
9.8 ± 2.9 



2.63 ± 0.01 2.63 fixed 
< 27 27 fixed 

1.1 ±0.3 < 1.8 
S xvi Lya, SixivRRC 



2.97 
< 29 
2.8 ±0.4 
ArKa 



0.01 



2.96 
< 38 
4.7 ±0.9 



0.01 



3.14 fixed 3.14 ±0.02 

56 fixed 56 ± 32 

< 0.69 4.3 ± 1.4 
ArxvuHea 

3.71 ±0.01 3.75 ±0.02 

< 28 < 46 
2.7 ±0.5 3.4 ±1.1 
CaKa 

6.40 ± 0.01 6.42 ± 0.01 

< 15 22 ± 10 
62.5 ±1.9 150 ±7 
Fe Ka I - xvn 



6.53 ± 0.04 
206 ± 30 
19.2 ±2.9 



6.72 ± 
< 42 
20 ±5 



0.02 



2.46 fixed 
20 fixed 

< 6.1 

2.63 fixed 
27 fixed 

< 2.0 

2.97 fixed 

29 fixed 

< 4.6 

3.14 fixed 
56 fixed 

< 0.81 

3.71 fixed 

30 fixed 

< 6.3 



+0.05 
-0.04 



6.41 
< 65 

5 3+ 73 

-2.1 



6.72 fixed 
42 fixed 
<0.92 



Fe Ka xvm - xxiv, Fe xxv Hea 



0.01 



7.06 
< 28 

11.3 ± 1.1 
Fe KyS I - xvil 



7.11 ±0.01 7.06 fixed 
37 ± 25 28 fixed 

60 ±7 < 13.1 



3.2.1 . The Fe-line complex 

Figs.|2]and|4]cleariy show that the Fe Ka line around 6.4 keV 
is the strongest line in the eclipse and egress spectrum (6.3 X 
1CT 4 and 1.5 x 1CT 3 photons cm" 2 s _1 , respectively). A quick 
look at the piled-up part of the lightcurve shows that the line 
is clearly detected throughout the orbit. However, in the low- 
flux spectrum the line has a strongly reduced flux (5 x 10~ 5 
photons cm" 2 s" 1 ; Fig. [5). In the eclipse and egress spectrum 
the measured line energy is consistent with an ionisation degree 
ranging from near-neutral iron to about Fe xvn. 

A second, broad line is detected in the eclipse spectrum at 
a slightly higher energy around 6.5 keV. Its energy is consis- 



tent with Fe Ka produced by more highly ionised iron (up to 
about Fexxiv). Since this line is very broad (cr ~ 0.3 keV) and 
close to the 6.4 keV line, we consi der a second possibility (as 
reported bv lLa Barbera et all200ll on BeppoSAX observations 
of LMC X-4) in which two gaussians at 6.4 keV are used to 
mimic one line with a narrow core and broad wings. This does 
not result in a good fit. In egress we do not find evidence for a 
second, broad line around 6.5 keV, but a narrow line at 6.7 keV 
is present, which is consistent with Fe Ka from almost com- 
pletely ionised, He-like iron (Fexxv). In the low-flux spectrum 
we do not detect this line. With its energy and width fixed to the 
values found for the line in the egress spectrum, we can place 
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Fig. 5. XMM/EPIC observation of the 4U 1700-37 during low flux. The soft excess has been fit by a third power-law component. 
Only some of the many lines observed in eclipse and egress are detected. Especially the strong reduction of the iron complex is 
remarkable. The RGS spectrum is displayed in the left bottom corner of the upper panel. The EPIC/MOS-2 residuals are shown 
in the middle panel, and the RGS residuals in the lower panel. 



an upper limit on its flux of 0.92 x 10~ 5 photons crrT 2 s _1 . This 
is much less than the measured flux in egress, similar to the 
strongly reduced flux of the Fe Ka line around 6.4 keV. 

In addition to the Fe Ka lines we detect the Fe K/3 line 
at 7.1 keV in the eclipse and egress spectrum. In the low-flux 
spectrum again only an upper limit can be placed on the flux 
of this line. Similar to the second Ka line, the energy of the 
Fe K/3 line is slightly higher in egress than in eclipse (7.11 and 
7.06 keV, respectively). This is consistent with iron ionised up 
to about Fe xvn, which is similar to the iron Ka line at 6.4 keV. 

The fluxes of the different iron lines show a similar be- 
haviour. All the lines are much weaker in the low-flux spec- 
trum than in the egress spectrum, and the fluxes of the near- 
neutral iron Ka and K/3 lines are higher in egress than in 
eclipse. However, there are some differences which become 
apparent when calculating the ratio Fe Kfi/Ka. In eclipse we 
find a value of K/3/Ka = 0.12*ggJ, which is in agreement 
with th e value of 0.13 that is expected for an optically thin 
plasma (Kaas tra"& Mewdll993l) . However, in egress this ratio 
increases to K/3/Ka — 0.35^"^. Due to the high upper limit on 
the flux of the K/3 line in the low-flux spectrum this ratio could 
be as high as 4.1. 



The presence of a K-edge from near-neutral iron at 
7.117 keV alters the shape of the continuum around this en- 
ergy. If the edge is not included in the model, this would result 
in an increased flux of the Fe K/3 line as well as the value of the 
K/3/Ka ratio. A deeper edge will thus decrease the K/3/Ka ra- 
tio. An edge is clearly present in both the eclipse and the egress 
spectrum (see Figs. |2 and 0}, which is taken into account by 
the absorption component of the model. We do not find evi- 
dence for a deeper edge. If we put the abundance of Fe at zero 
the edge is not modelled by the absorption component. If we 
add an edge to the model separately, this results in an edge at 
slightly higher energy (7.17+° ^ keV). The flux of the K/3 line 
reduces to ~ 3 x 10~ 4 photons crrT 2 s _1 . The K/3/Ka ratio then 
beco mes ~ 0.15, which is c onsistent (within its errors) with 
0.13 dKaastra & Mewell993h . 



3.2.2. Line flux variation with orbital phase 

We compare the flux from eclipse to egress of the different dis- 
crete recombination lines considering only the lines that we can 
clearly separate from the other lines and for which we have a 
unique identification. The fluxes of Si xiv Lya, Si xm Hea and 
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Fig. 6. The flux ratio of the lines in eclipse and egress as a func- 
tion of the ionisation parameter log Lines corresponding to 
higher £ values become stronger going from eclipse to egress. 



Ne x Lya remain constant within their errors, while the fluxes 
of S xv Hea, Ar xvn Hea and Fe xxv Hea increase towards 
egress (see Fig.|6). 

We use a power-law model with phot on index a — 1 .4 in 
the XSTAR photoionisation code jBautista & Kallmanl feoOll) 
to determine the ionisation parameter £ (see Sect. |4ji for the 
reported ions, assuming a photoionised plasma. The ionisation 
parameter £ is defined as 



(1) 



where Lx is the bolometric X-ray luminosity, n the number 
density at the dist ance r from the i onising source in an opti- 
cally thin plasma (Tarte ret all 19691) . The ionisation parameter 
then determines the local degree of ionisation of the plasma 
and can be used to predict the population of ions expected for 
a certain value of this parameter. 

The values of £ are listed in Tabled We apply these values 
on the selection of lines and plot the ratio of the line fluxes 
in eclipse and egress as a function of £ (see Fig. |5J. We find 
that the lines with log £ > 2.7 are stronger in egress while the 
lines with a lower value for the ionisation parameter remain 



Table 4. The value of the ionisation parameter £ corresponding 
to the reported emission lines. 



H-like 


He-like 


ion 


log£ 


ion 


logf 


vm 


1.6 


Ovn 


1.0 


Nex 


2.2 


Neix 


1.6 


Mgxn 


2.6 


Mg xi 


2.1 


Si xiv 


2.8 


Sixm 


2.4 


S xvi 


3.2 


Sxv 


2.6 


Arxvm 


3.4 


Arxvn 


2.8 


Fe xxvi 


4.2 


Fexxv 


3.6 



constant, as expected for a plasma that is photoionised by the 
X-rays from a central source. 



4. Discussion 

We interpret our observations including information obtained 
in other wavelength domains. The general picture we have in 
mind is that of a compact object, which accretes from the ex- 
tended stellar wind of the O supergiant companion in which it 
is embedded. Due to this accretion an X-ray continuum is pro- 
duced, which photoionises its near surroundings (see Fig. 0. 
The degree of ionisation (expressed in terms of in this re- 
gion decreases with distance from the X-ray source. Emission 
lines corresponding to higher £ values originate closer to the 
X-ray source (see. Eq.Q. 

4.1. Continuum emission 

We have shown that the high energy part of the continuum 
of 4U 1700-37 can in all cases be represented by a hard, di- 
rect component and a (Thomson) scattered component, which 
are modelled as two power laws with identical photon index, 
but different absorption coefficients. The amount of absorption 
on both components depends on the line of sight to the X-ray 
source and thus on the orbital phase. 

The direct and scattered component can be clearly sepa- 
rated in the egress spectrum, when the absorption of the di- 
rect component is high. The unabsorbed flux (Fx) in the energy 
range 0.3-1 1.5 keV of the direct compone nt is ~ 5 x 10~ 9 er g 
s _1 cm -2 . Adopting a distance of 1.9 kpc jAnkav et alJ EoOl ). 
this yields an unabsorbed luminosity Lx ~ 2.2 x 10 36 erg s _1 , 
consistent with previously reported values. 

During eclipse there is still evidence for the presence of 
the direct component, although its flux is strongly reduced 
(F x ~ 1 X 10~ 10 erg s -1 cirT 2 , which gives L x ~ 5 x 10 34 



erg s '). Note that our dataset does not cover a full eclip se, but 
starts a t cj> ~ 0.08 accordi ng to the ephemeris of Rubin et all 
dl996l) . lRubinetalJ(ll996l) derive an eclipse semi-angle 8 e = 
28°. 6 + 2°.l, using BATSE on Compton GRO, which has an 
energy range of 20-120 keV. Therefore, the "physical" eclipse 
should end already at <p ~ 0.08. However, at soft energies the 
eclipse lasts longe r due to the strong a bsorption of soft X-rays 
by the stellar wind. lHaberl et alJ i 19891) found e = 44°.6 + 0°.6 
in EXOSAT observations in the 2-10 keV energy range. This 
results in an "end of the eclipse" at ~ 0.12, consistent with 
our observations. Due to this energy dependent definition of 
eclipse, a small contribution to the continuum with high ab- 
sorption can still be expected. 

We modelled the low-flux spectrum at orbital phase (f> ~ 
0.25 in two ways, differing in how the high-energy part of 
the spectrum is fitted. In model A (see Table |2) a highly ab- 
sorbed direct component and a scattered component are used, 
with similar unabsorbed fluxes (Fx ~ 5 x 10~ 10 erg s cm -2 ; 
Lx ~ 2 x 10 35 erg s ). Model B (see Table[2} contains only a 
direct component with an unabsorbed flux similar to that de- 
rived for model A. This component is modelled as a power 
law modified by a high-energy cutoff and is the same model 
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Fig. 7. This sketch of the geometry of 4U 1700-37 includes the 
size of the region where q < 4, the upper limit on t he size of 
the ionised region reported by Ivan Loon et al.l J200ll) 



iRevnolds et alJ i 19991) used to fit their BeppoSAX observa- 
tions. For both models the flux of the direct component is a 
factor ~ 9 lower than in egress. In model A the absorption on 
the direct component is relatively high (Nh ~ 2.0 x 10 23 cm 2 ), 
whereas in model B the absorption is lower (Nh ~ 7 x 10 22 
cm" 2 ). This value i s consistent with calculations reported by 
IHaberl et all dl989h . who predict N H ~ 5 - 7 x 10 22 cirT 2 at 
orbital phase <p ~ 0.25. 

If we compare the flux of the direct component in egress 
to that in low flux, we find that a similar decrease in the scat- 
tered component should result in a non-detection. A strong re- 
duction of the flux of the direct component and consequently 
of the scattered component as inferred from model B can be 
explained by a reduction in accretion, suc h as the case in a 
"clumped" wind structure JSako et all (f 1999) also find evidence 
for clumping in the stellar wind of the optical companion of 
Vela X— 1. In case of Bondi-Hoyle accretion, often used to ex- 
plain th e prod uction of X-rays in wind-accreting HMXBs (e.g 
Nagase 1989), the amount of material accreted (and thus Lx) 
scales as p/v 3 , where p is the density of the accreted material 
and v its velocity. Therefore, the accretion is very sensitive to 
changes in the density- and (especially) in the velocity structure 
of the stellar wind (see Kape r et alJl993l) . The X-ray lightcurve 
(see Fig.Q shows strong variability, which suggests a contin- 
uously changing accretion rate onto the compact object. Given 
the duration of the low-flux interval (~ 5 ks) and assuming a 
typical wind velocity v ~ 1000 km s" 1 the region shou ld ex- 
tend upto ~ 5 x 10" cm (i.e. ~ 0.4 R+: IClark et all2002h . 

Model A can only be consistent with this scenario if the 
scattered component is not affected by the drop in accretion, 
since this model predicts a scattered component with a similar 
unabsorbed flux as the direct component. This is only possible 
if it is produced sufficiently far out in the system, that it does 



not notice a change in accretion rate for ~ 5 ks. This would 
require a distance of ~ 10 14 cm (~ 70R+), whereas the orbital 
separatio n of the compact object and the optical companion is 
~ 2 R* iHean & Corco'rIr] IT992). This is not consistent with 
the relatively high column density (Nh ~ 7 x 10 22 cm" 2 ) toward 
the scattered component. 

A possible alternative scenario that can be described by 
model A is the presence of an accretion wake, formed in front 
of the compact object, which creates a kind of eclipse that 
blockes the photons originating from the direct surroundings 
of the compac t object. The prese nce of such a wake has been 
suggested by Hab erl et alJ i 19941) . based on a difference in the 
temperature of the bremsstrahlung component found in Ginga 
observations. In this scenario the interval of low flux should 
be locked to orbital phase, which is not consistent with other 
observations that report on low-fl ux intervals at different or- 
bital phases. Boros on et al find a quiescence period in 
Chandra observations around <p ~ 0.68. Therefore, model A 
does not give a good physical description of the low-flux spec- 
trum. 

In the case of an exponential cutoff, the values of E c and 
Ef can not be determined very well, due to the energy range of 
the EPIC/MOS-2 camera, but using fixed values for E c and Ef 
in eclipse and egress, model B can be used to describe all three 
spectra. There is an indication for an increase of the photon 
index from eclipse and egress to the low-flux interval. Although 
the value of the photon index is correlated to the value of Ef, 
one can not exclude a change in the physical conditions of the 
accretion during the low-flux interval. 

We conclude that model B probably best describes the 
physical conditions of the low-flux interval, in which the scat- 
tered component is not detected. A broader energy range is 
needed to further constrain the parameters of the high-energy 
cutoff. 

4.2. Soft excess 

A soft excess is clearly present in all spectra and can be mod- 
elled with components of different physical nature, modified 
by an absorption which is consistent with the interstellar value 
for this so urce. The previously reported bremss trahlung com- 

I '— | 1 r- 1 <- II 1 ° 

ponent (Haberl & Day 1992; Haberl et al. 1994) results in a fit 
with a temperature that can not be constrained and a shape that 
mimics a power law. Also in the analysis of previously ob- 
tained data, the tem perature of the bremsstrahlung componen t 
had to be fixed (e.g. lHaberl et al.ll 19941 IRevnolds et alJll999T) . 
Therefore, we do not obtain more information on the physical 
origin of this component than with a simple power law model. 

Another possible model for the soft excess is black- 
body emission. This component is often visible in HMXBs 
hosting an accretion disc (e.g. LMC X— 4 a nd Cen X— 3; 
lLa Barbera et alJEoOll and iBurderi et alJl2000l respectively). 
We can calculate the size of the emitting surface of the black- 
body component: 

R km = 3.05 x 10Vp C F - 5 r ke 2 v , (2) 

in which R km is the radius of the emitting surface in km, d\^ c 
the distance to the source in kpc, F the flux in erg s" 1 cm" 2 and 
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7\ e y the t emperature in keV. Adopting a distance of 1.9 kpc 
(Ank avetalihoOll) . we find an upper limit for the radius of 
the emitting volume of 52 km for the low-flux spectrum. For 
the eclipse and the egress spectrum we find upper limits for the 
radius of 0.13 and 0.34 km, respectively. These latter values 
are very small. Adding the fact that the absorption is very low, 
the blackbody component can not originate from an accretion 
disc or (close to) the compact object. We therefore conclude 
that although the soft excess can be described by a blackbody 
component, a proper physical interpretation is lacking. 

The last possible continuum model which fits the soft 
excess is scattering. We examined both dust scattering and 
Thomson scattering. Dust scattering is energy dependent (oc 
1/Zs 2 ), which results in a steeper power law, mainly influenc- 
ing the high-energy part of the spectrum. The difference be- 
tween dust- and Thomson scattering can only be measured if 
the high energy part of the spectrum is not dominated by the 
strong direct component, i.e. during eclipse. Using dust scat- 
tering in eclipse does not improve the fit, whereas Thomson 
scattering does. Furthermore, this component has to originate 
from the outskirts of the system due to its low absorption. This 
is not consistent with a continuum that is absorbed in the ex- 
tended stellar wind on its way out, before it is partly scattered. 
Therefore, we conclude that although Thomson scattering is a 
model that fits the data well in all cases, it is unlikely to be the 
physical origin of the soft excess. 

Recently, iBoroson et alJ d2003l) reported on Chandra ob- 
servations of 4U 1700-37. They claim that the soft excess can 
be explained by a blend of lines only. However, the statistics of 
their observations are low at these energies and they only report 
one unresolved line below 1 keV. We tried to fit the soft excess 
in our spectra in a similar way by removing the soft contin- 
uum component and using broad gaussians to create a pseudo- 
continuum. Although this results in stable fits, the lines are 
4-5 times broader than the energy resolution of EPIC/MOS- 
2 at low energies. Furthermore, the absorption of the direct 
and scattered component is ~ 3 times lower in eclipse than in 
egress. We conclude that this model does not provide a satisfy- 
ing physical description of the continuum. However, it does not 
exclude the possible presence of a blend of lines. Therefore, it 
remains a matter of debate whether the soft excess is real or 
due to a blend of unresolved emission lines. 

4.3. The Fe-line complex 

The line spectrum consists of fluorescence emission lines of 
near-neutral species and discrete recombination lines of He- 
and H-like species. In eclipse and egress many lines are de- 
tected, while the low-flux spectrum shows only a few. Many 
lines are blended with each other or have multiple identifica- 
tions, which complicates the analysis. 

One of the most striking results is that the iron Ka com- 
plex around 6.4 keV is strongly reduced in flux in the low-flux 
spectrum (see Fig. [5}, while it is clearly visible in the rest of 
the (piled up) lightcurve. This is consistent with the reduced 
flux of the continuum above 6.4 keV and therefore another ar- 



gument in favor of the scenario invoking a reduced accretion 
rate during the low-flux interval. 

Moving from eclipse into egress, we find an increase by 
a factor ~ 2 of the flux of the Fe Ka line at 6.4 keV. This 
suggests that the formation region of the fluorescence lines 
from near-neutral elements is not far from the compact object. 
However the presence of discrete recombination lines from He- 
and H-like elements indicates a highly ionised region. In some 
systems an accretion disc may be the source of t hese emis- 
sion lines (e.g. LMC X-4: La Barbera et al1l20Q lh . However, 
4U 1700-37 is a wind accreting system, for which no evidence 
for an accretion disc has been found. Another possible scenario 
(e.g. ISako etalj Eool) relates to the clumping of the stellar 
wind. If the wind is inhomogeneous, cool dense regions can 
form the near-neutral material, responsible for the presence of 
the fluorescence lines. This would result in a hybrid medium of 
cool dense and highly ionised regions. 

The expe cted KJ3/Ka ratio for an optically thin plasma 
jKaastra &Mewell 1993b should be around 0.13. In Sect. EI] 
we show this is the case for eclipse, but in egress this is 
only true if the Fe edge is shifted to a slightly higher energy 
(7.17+0Q2 keV). An edge at such ener gy can be p roduced by 
slightly more ionised iron, upto Fe m l lLotzll 19681) . Ionisation 
stages of iron like Fern and Feiv are abundant in the stel- 
lar winds of OB supergiants and are important wind drivers. 
Although we can exclude the edge to be produce d by Feiv 
or hig her, this is not consistent with the results of IVink et alJ 
(2000) who find that the stellar winds of OB-stars with tem- 
peratures higher than 25000 K are line driven by Fe iv, which 
recombines to Fe m at larger distances from the stellar p hoto- 
sphere, typically around R ~ 6 - 8 R* ( Vink et al. 1999). This 
is bey ond the orbital separation of ~ 2 R* ( H eap & Corc oran 
Il992l) and thus outside the region probed by the compact object 
during the egress interval (see Fig. [7}. 

4.4. Recombination lines 

We detect discrete recombination lines in both the eclipse and 
the egress spectra. The emission lines reported in Sect. l3.2l can 
be identified with He- and H-like ions, with ionisation param- 
eters in the range 2 < log £ < 4. If we compare the flux of 
the lines that are clearly separated and have a unique identifi- 
cation, as a function of log we see a systematic change in 
flux of some of the lines from eclipse to egress. The elements 
with log £ > 2.7, i.e. S xv Hea, Ar xvn Heo- and Fe xxv Heo-, 
grow in strength (see Fig.[6jl, while Si xiv Lye, Si xin Hea and 
Ne x Lya remain constant within their errors. The observation 
of discrete recombination lines in an eclipse spectrum directly 
implies that the formation region extends beyond the size of the 
O supergiant. In our observations the eclipse interval covers a 
range in orbital phase of 0.07 < <f> < 0.14; therefore, the size 
of the ionised region should be ~ 0.5 R* (see Fig. [7}. Because 
the lines with log £ > 2.7 are stronger in egress, the Stwmgren 
zone, often defined as an ionised region with log £ ~ 3, can not 
be too large. 

At these relatively high values of £, inside this ionised 
region ions like Civ, Siiv and Nv should be totally ab- 
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sent, resulting in an observable modulation of the UV res- 
onance lines formed by these ions. This orbital modulation 
of UV resonance lines has be en predicted for 4U 1700-37 
bv lHatch ett & M cCravl i 1 9771). but has never been observed 
(Du pree et al.l Il978l iKaper et al.lll993l) . A possible explana- 
tion for the absence of the Hatchett-McCray (HM) effect in 
4U 1700-37 is that the stellar wind has a non-monotonic 
("turbulent") velocity structure which can hide the presence 
of an ionisation zo ne from detection in an UV resonance line 

jKaneret alii 9931 

Ivan Loon et alJ(l200lh derive the size of the ionisation zone 
for a few HMXBs, based on a wind model taking the non- 
monotonic velocity structure into a ccount. The size of the ion- 
isation zone can also be expressed (Hatchett & McCrav 1977) 
in terms of the non-dimensional variable q, that is defined as 



q = 



{n x D 2 



(3) 



where ^ is the ionization parameter, Lx the X-ray luminosity, 
rtx the number density at the orbit of the ionising source and D 
the distance from the center of the O supergiant to the ionising 
source. Substitution of eq.^in eq.[3]then yields 



n x D 2 



(4) 



In a medium with constant density, each value of q corresponds 
to a circle centered on the X-ray source, where the radius in- 
creases with decreasing value of q. In a stellar wind, the den- 
sity falls off roughly as r+ 2 , such that surfaces of constant q 
are small and closed (high q value; see Fig. 01 or larg e and 
open with the ionising source off-center (low a: Ivan Loon et alJ 
120011) . To each value of q corresponds a value of i.e. the 
higher the d egree of ionisation, the higher the value of q. For 
4U 1700-37 lvan Loon et aD J200ll) set an upper limit of q ~ 4 
to the size of the ionisation zone based on the absence of the 
HM effect in HD 153919 (see Fig.Q. They predict that q must 
be of the order 200, for log £ = 3. Such a region would be much 
smaller than the size of the O supergiant. The detection of dis- 
crete recombination lines with high £ values during the eclipse 
spectrum indicates that the size of the ionisation zone must be 
larger, which is still consistent with the reported value for the 
upper limit, q < 4 Jvan Loon et alJ2o"oil) . 



5. Summary and conclusions 

XMM-Newton observations of 4U 1700-37 during eclipse 
(0.07 < 4> < 0.14) and egress (0.14 < (f> < 0.17), as well as 
a low-flux interval around (f> ~ 0.25, produced the following 
results: 

- The continuum is well represented by a combination of two 
power-law components: (i) a direct component produced 
by the X-ray source and (ii) a (Thomson) scattered compo- 
nent originating in an extended region surrounding the X- 
ray source. The extent of this region must exceed the size 
of the O supergiant as the scattered component is detected 
during eclipse. 



- The low-flux interval is likely due to a lack of accretion 
onto the compact object due to a varying density and/or ve- 
locity structure of the accreting material, such as expected 
for a clumpy stellar wind of the optical companion. 

- At low energies (< 2 keV) a soft excess is present, similar 
to what has been found in earlier observations with, e.g. 
EXOSAT and Ginga. We can exclude a thermal origin in 
the form of bremsstrahlung or blackbody emission. 

- In all spectra we detect fluorescence line emission from 
near-neutral iron around 6.4 keV, probably produced by 
dense "clumps" in the stellar wind close to the compact 
object. In the low-flux interval this line is strongly reduced 
in flux in tandem with the continuum above 6.4 keV. 

- The initially measured high Fe KJ3/Ka ratio in egress can 
be explained by an Fe edge due to moderately ionised iron, 
upto Fe m. This is consistent with the ionisation degree of 
iron further out in the stellar wind of the optical companion, 
but not with the region probed by the compact object in 
egress. 

- We detect many discrete recombination lines from He- and 
H-like species in both eclipse and egress. Some of the lines 
can be isolated, while most of them are blended with each 
other. The presence of these lines in eclipse directly implies 
that the formation region has a minimum size of around 
0.5R*. Some of the lines with a relatively high value of the 
ionisation parameter ^ grow stronger towards egress. This 
implies that the ionisation zone can not be very large either. 
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